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Abstract 

We consider stationary versions of two discrete variants of Hammersley’s process in a 
finite box. This allows us to recover in a unified and simple way the laws of large numbers 
proved by T. Seppalainen for two generalized Ulam’s problems. As a by-product we obtain 
an elementary solution for the original Ulam problem. 

We also prove that for the first process defined on Z, Bernoulli product measures are the 
only extremal and translation-invariant stationary measures. 

MSC 2010 Classification: 60K35, 60F15. 
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1 Introduction 

In a celebrated paper, J.M.Hammersley used Poissonization to attack the so-called Ulam problem 
of the typical length i{n) of the longest increasing subsequence of a uniform permutation of size 
n. Namely, he reduced this problem to finding the greatest number of points of a Poisson point 
process inside a square, an increasing path can go through. He proved (|3j. Theorem 4) that 
i{n)/\fn converges in probability to some constant c, sometimes refered to as the Ulam constant, 
and conjectured that c = 2. 

The proof of c = 2 was achieved independently by Logan and Shepp and by Vershik and 
Kerov in 1977, using algebraic methods. Various authors were then interested in finding a more 
probabilistic proof of this result. First, Aldous and Diaconis [T] gave one, using the properties 
of what they called Hammersley’s process, which was implicitly introduced in [1] (p.358 and 
following). Hammersley’s process is continuous in time and space and Aldous and Diaconis 
studied its properties on the infinite line, in particular its stationary distributions. A few years 
later, Groeneboom [3] and Gator and Groeneboom [2] studied Hammersley’s process on a quarter 
plane. By adding what they called Poisson sinks and sources on the x and y-axis, they also 
found a stationary version of this process on the quarter plane. Using this point of view, they 
were able to recover again the value of c. 

In this paper, we study two discrete variants of Ulam’s problem. Namely, for all p in [0,1], 
n, m > 1 , we replace the original Poisson point process by the following random set ^ of integer 
points: each integer point of the rectangle [1, n] x [1, m] is chosen independently with probability 
p. We are interested in the two following quantities. 
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The length of the longest increasing subsequence through points of 


= max{L; there exists • • •, (u, ji) G < • • • < and ji < j2 < • • • < Jl} ■ 

• The length of the longest non-decreasing subsequence through points of 
(‘ 2 '] 

^\n,m) = max{L; there exists • • •, (u, ji) G < • • • < and ji < j2 < • • • < ji} ■ 


One aim of the present work is to recover by simple and unified probabilistic arguments 
the first order asymptotics of and already obtained by T. Seppalainen in two 

independent papers, for in [9] and for in [10]. In the following result the variables 

. ) are coupled in the obvious way. Moreover, for any a, 6 > 0, . n stands for 

V (an,om) 

T (0 

^ilan\,lbn\)- 

Theorem 1. For a,b > 0, we have, when n tends to infinity 


(an,bn) 

n 

^( 2 ) 

(an,bn) 

n 


y/p{2^/afi - (g -F b)y/p) 

< 1 — p 

min{a, b} 

K 

j2^abp(l - p) + {a- b)p 

I a 


if p < min{a/ 6 , b/a}, 

( 1 ) 

otherwise, 

ifp < a/{a + b), 
otherwise. 

( 2 ) 


To prove this, Seppalainen associates to each problem a particle system on the infinite line 
Z and checks that Bernoulli product measures are stationary. The position of particles at a 
given time is then characterized as the solution of a discrete optimization problem, and the a.s. 
limit of identified using convex analysis arguments. Note also that for the second 

1—1 ( 2 ^ 

question, Johansson used later ([ 6 |, Th.5.3) another description of Lf^Jnbn)/^ proved that 
the rescaled fluctuations converge to the Tracy-Widom distribution. 

In | 8 |, Rolla-Sidoravicius-Surgailis-Vares considered the closely related model of last-passage 
percolation with geometric weights. Here is a quick description of the model. Let {^{i, j))i<i<n,i<j<m. 
be a family of independent geometric random variables with parameter p. Set: 

=™ax|^e(7(i))| 

where the maximum is taken over all paths 7 from (1,1) to (n, m) with (1,0) or (0,1) steps. With 
respect to or we have simply changed the weight of each point (from Bernoulli to 

geometric) and the set of paths considered. Rolla-Sidoravicius-Surgailis-Vares gave a new proof 
of the following result: 

r(3) 

(n,n) a.s. ^ ‘^y/P 

n n ^+00 1 — y/p 

(The original proof was given by Jockusch-Propp-Shor in |5|.) 

Our proof of Theorem 1 is inspired by the proof of Rolla-Sidoravicius-Surgailis-Vares of the 
above result. The strategy is to consider a particle system on a bounded domain which turns 
out to coincide with the restriction of Seppalainen’s particle system on Z. This simplihes the 
definition of the process (especially, in the case of the second problem). Moreover, it turns 
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out that a local balance property around a single site, see Lemma and below, is enough to 
check the stationarity of the process. Our proofs are essentially the same for both models. This 
kind of remarkable local balance property also occnrs in last-passage percolation with geometric 
weights (see Rolla-Sidoravicins-Surgailis-Vares , Section 3.1 in [8j or Seppalainen, Lemma 2.3 
in m)- The particle systems studied here provide two other examples where such a property 
holds. Theorem 1 is then proven by investigating the behaviour of the models under the different 
stationary measures which have been exhibited. 

A nice by-product of our proof is that we obtain non-asymptotic estimates which provide 
an elementary proof of c = 2 for the original Ulam problem (see the discussion at the end of 
Section . 

In the last section, which can be read independently, we complete the stndy of the original 
infinite discrete Hammersley process introdnced by Seppalainen [Hj for the Problem 1. We prove 
that the only extremal translation-invariant stationary measures of this process are the Bernoulli 
product measures. 

2 Discrete Hammersley’s processes 

Like Hammersley did, we construct two sets of broken lines whose cardinality is respectively the 
variable n and For each case, we first introduce a partial order: 

(n,mj (n,m) ’ ^ 

Case 1: Case 2: 

(x,y)-< (x',?/') iff {x,y) < iS 

X < x' and y < y' x < x' and y < y' 

Now, Hammersley’s lines are paths starting from the top side of the rectangle [l,n] x [l,m], 
ending at its right side and making only south or east steps. They are constructed recursively. 
The first line is the highest non-increasing path connecting the minimal points of ^ for We 
withdraw these points from ^ and connect the new minima to get the second line, and so on. 
In the below picture, n = m = 8 and = 4, = 5 (crosses denote points of 

Hammersley’s lines are in bine, one of the longest subsequences is shown in a red dashed line): 




Problem 1 Problem 2 

We claim that the number H; ' , of Hammersley’s lines is equal to L; ’ Indeed, note that, 

(n,m) ^ ’ ’ 

(i) 

in both cases, is equal to the length of the longest increasing subsequence trongh points 

of ^ for the partial order Each Hammersley’s line connecting minimal points of a subset of 
the points of a maximizing seqnence necessarily lie on distinct Hammersley’s lines and 

(i) (i) 

therefore To prove the opposite inequality, we observe that Hammersley’s lines 
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give a construction of an increasing sequence of points ai,... a (i) of ^ in the following way. 

(n,m) 

For a (i) we take any point of ^ which lies in the top right Hammersley’s line. Then, for each 

i, ai-i is taken as one of the points of ^ which lies in the I — 1-th line and such that a£_i -< a^. 
The existence of such a^-i is granted. Indeed, if no such point existed, would have belonged 
to the {I — l)-th line. This proves w 

Let us note that the construction of lines above is consistent in n and m, i.e. for n' < n 
and m! < m the restriction of Hammersley’s lines to the rectangle \^,n'\ x [1, m'] coincides with 
Hammersley’s lines constructed from the points of ^ in this smaller rectangle. In particular, for 
n > 1, we can define the two discrete time processes {X^)t>o := x G [[l,n]])f>o and 

{Xi)t>o ■= {Xf{x), X G [l,n]])i>o defined by 


^t(^) = 


1 if there is a Hammersley’s line on the vertical edge {(x, t), (x, f + 1)} 
0 otherwise. 


( 3 ) 


Note that on each vertical edge there is at most one Hammersley’s line. Moreover, in Model 1, 
there is also at most one Hammersley’s line on any horizontal edge and Hammersley’s lines do 
not intersect each others. We will say that there is a particle at time t at position x in Model 
i if Xl{x) = 1. Hence, for both models we start with the initial condition Xq = 0 and one can 
check that both processes are Markovian. Besides, since Hammersley’s lines start all from the 
top side of the rectangle, 

for both models, the number of particles at time t is equal to 

In both models, the dynamic of the particle system is quite simple. As the description of 
the dynamic is not needed in the proof of Thereom we do not give such a description here. 
However, we refer the interested reader to Appendix [A} 


3 Sinks and sources 


Since, in both models, the number of particles can only increase as time goes on, both processes 
converge to their unique stationary measure := 1. In this section, we are going to modify 
a little bit the problems we are looking at, such that the particle systems associated to the 
new problems admit less trivial stationary measures. To do so, we must first introduce some 
notation. 

From now on, for any integer x > 1 and t > 1, we define the random variable ^t(x) by 


it{x) 


1 if (x, f) G ^ 
0 otherwise. 


ii) 

With this notation, we can write the quantities as 


L 


(i) ^ 

{n,m) 


max I 


( 4 ) 


( 5 ) 


where the maximum is taken over all the discrete paths (x(-), t(-)) in [1, n] x [1, m] starting from 
(1,1) which make only steps in N x N for Model 1 and in N x Z_|_ for Model 2. Given two 
sequences of integers and (^o(3;))a:e[i,n] j we now consider the quantities 


{n,m) 



( 6 ) 
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where the maximum is taken over all the discrete paths (x(-), t(-)) in [ 0 , n] x [ 0 , m] starting from 
( 0 , 0 ) such that: 

• There exists an integer u G Z_|_ such that: either they first make u (0,1) steps or they first 
make u ( 1 , 0 ) steps. 


• Then, they only make steps in N x N for Model 1 and in N x Z+ for Model 2. 

Using the same terminology as in [2], we say that there is k sinks (resp. sources) at site (0, t) 

(resp. (x,0)) if ^t(O) = k (resp. ^o{x) = k). The sources will be drawn on the x-axis and the 

(i) 

sinks on the y-axis. Hence, in other words, is equal to the maximum number of points 

an increasing (resp. non-decreasing) path can go trough when the path can pick either sources 
or either sinks before going into N x N. The law of the sources and sinks will be specified later, 
let us just say that the sources will take values in { 0 , !}”■ whereas the sinks will take values in 
{0,1}"* for Model 1 and in Z™ for Model 2. 




Model 1, , = 5 Model 2, . = 7 

’ (71,nj ’ (77,7lJ 


ii) 

As before, we can construct a set of broken lines whose cardinality is equal to Now, 

Hammersley’s lines are paths starting from the top side or the left side of the rectangle [0, n] x 
[0, m], ending at its bottom side or its right side and making only south or east steps. They are 
still constructed recursively. Assume that the sink of minimal height is located at ( 0 , ti) and 
the leftmost source at (xi,0). Then, the first line is the highest non-increasing path starting 
from ( 0 ,ti), ending at (xi, 0 ) and connecting minimal points {x,t) of ^ for the partial order ^ 
introduced in the previous section such that (xi,0) {x,t) and (0, ti) 7 ^ {x,t). If there is no 

sink (resp. source) at all, we do the same construction except that the line starts from the top 
of the box (resp. ends on the right of the box). Then, we withdraw these points from one of 
the sink located at ( 0 ,ti) and the source located at (xi, 0 ) and we do the same procedure with 
this new set of sources, sinks, and set of points in [[l,n]] x to obtain the second line. The 

algorithm goes on until there is no more points, sinks and sources. Using the same argument 
as in the previous section, one can easily prove that the number of lines obtained with this 
procedure is equal to 

With a slight abuse of notation, we still denote {Xl)t>o := {Xl{x), x G [[l,n]])j>o and 
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Xlix) = I 

Hence, the 


(X|(x), X G [[l,n]])f>o the two discrete time system of particles defined by 

1 if there is a Hammersley line (with sources and sinks) on the edge {(a:, t), (x, t + 1)} 
0 otherwise. 

processes now start from the configuration given by the sources i.e Xq{x) = io{x). 


Theorem 2 (Stationarity of Hammersley’s processes on a bounded interval, with sources and 
sinks). For all n, 

• Model 1. For allp,a G (0,1), the process {X^{x),x G [[l,n]])i>o is stationary if sources 
are i.i.d. Ber(a) and sinks are i.i.d. Ber(a*) (and sinks independent from sources) with 


* p(l - a) 

' ~ - 7 - 7 ' 

a p{l — a) 


( 7 ) 


• Model 2. For all p,a G (0,1) such that a > p, the process {Xf{x),x G |{l,n]])t>o is 
stationary if sources are i.i.d. Ber(a) and sinks are i.i.d. (Geo(a*) — 1) (i.e. P(^(0,t) = 
k) = Q!*(l — a*)^ for k = 0,1,...) (and sinks independent from sources) with 


a — p 
«(1 - P )' 


( 8 ) 


(Note that a* G (0,1) if a > p.) 

Remark 1. In f3i (resp. in JT^). Seppdldinen associates to Model 1 (resp. Model 2) a particle 
system defined on the whole line Z. One can check that the restrictions of these systems to [[1, nj 
have respectively the same transition probabilities that our processes X), i = 1,2 when sinks and 
sources are distributed as in Theorem Hence, 

• For Model 1 (resp. Model 2), Theorem^ should be compared to Lemma 2.1 in (resp. 
Proposition 1 in m) which states that Bernoulli product measures are stationary for the 
analogous model on the infinite line. 
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• The condition a > p in Model 2 reminds the condition on the density of the initial config¬ 
uration needed to define the process on Z in m- 

• The theorem states stationarity in time (^i.e. from bottom to top in our figures). In fact, 
the proof also shows stationarity from left to right. 

Let us collect a few consequences for further use. Let be the number of Hammersley’s 

lines in Model i with parameters a (sources), /3 (sinks) in the box [0,n] x [0,m] and let 
be the number of Hammersley’s lines that leave the same box from the top. Then, for i = 1, 2 
and every a, (3 G [ 0 , 1 ] 


ln,m) (n,m) {n,ra)'’ 

^(n mf ~ {sinks between 1 and m} . (9) 

This relies on the fact that any lines has to exit by the top or by a sink. Note that in the 

right-hand side of ([^ the two terms are not independent. Moreover, using the interpretation of 

^\nm) ^(n mf length of a longest subsequence {cf. ([^ and ([^), we have 


^(i) ^ ^(i),OL,l3 

{n,m) — (n^m) 


( 10 ) 


Besides, being by definition equal to the number of particles at time m, Theorem 

implies that 

Binom(n, a). ( 11 ) 

Proof. We aim to prove Theorem We first prove for both models a local balance property, in 
Lemmas 1 and 2. These lemmas are elementary but do not seem to be written elsewhere, and 
form the heart of our proof. 

Model 1. 

We first focus on what happens around a single point {x,t). Recall that ^tix) is a Bernoulli(p). 
Denote by X (resp. Y,X',Y') be the indicator that a Hammersley’s line hits {x,t) from the 
bottom (resp. left,top,right). 

Lemma 1 (Local balance for Model 1). The random variables {X',Y') are measurable with 
respect to {X,Y,f^tix)) and we have 


{X, Y, ftix)) = Ber(Q;) (g) Ber(Q;*) (g) Ber(p) 1 f {X', Y') = Ber(Q:) (g) Ber(a*) 


where a* is as in 0. 


X' ~Ber(a) 

^ ind. 

Y ~Ber(a*) . . Y' ~Ber(a* 

. St\T) 

\ ;* 
ind. 

^ X Ber(a) 
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Proof of Lemma^ Recalling that, in Model 1, Hammersley’s lines are non increasing lines which 
do not touch each other and noticing that X' = Y' = 1 X = Y = 0 and = 1, we get 
that 

r (X,y) ifx/y 

{X',Y') = l (1,1) if X = y = 0 and ^t(x) = 1. 

[ (0,0) otherwise. 

Hence, we see that {X', Y') G (t{X, Y, and one can easily check that if {X', Y') ~ Ber(a) (g) 

Ber(a*) with a* is as in (Q, then (X', y') ~ Ber(a) ® Ber(a*). □ 

We now explain how Lemma proves Theorem for Model 1. For every (x, 1) define the 
corresponding (X'(x),y'(x)) be the output of (x, 1). Independence of sources and sinks and 
Lemmaensure that the output of (1,1) is independent of sources at (2,0), (3,0),... and sinks 
at (0, 2), (0,3),.... Then a simple induction on x proves that the random variables (X'(x), x G 
[[l,n]]) are distributed as Ber(a). 

Model 2. The strategy is similar: 

Lemma 2 (Local balance for Model 2). The random variables (X',Y') are measurable with 
respect to (X, y, ^t(x)) and we have 

(X, y, ^t(x)) = Ber(a) (g) (Geo(a*) — 1) (g) Ber(p)^ ^(X', Y') = Ber(a) (g) (Geo(Q;*) — 1) 
where a* is as in Q. 

X' ~Ber(Q;) 

^ ind. 

y ~Geo(a*) - K , , Y' ~Geo(a*) - 1 

, ^t{x) 

\ / 

ind. 

^ X Ber(Q;) 


Proof of Lemma^ Recall that X, X' G {0,1} whereas y, y' G TP. Besides, using that Ham¬ 
mersley’s lines are non-increasing, the balance between incoming and outcoming lines at a site 
yields X' -|- y = X -|- y'. Moreover, X' = 1 iff (^t(x) = 1 or (X, y) = (1, 0)). 

We observe that the different cases can be summed up in 


(x',y') = 


( 1 , 0 ) 


if (x,y) = (i,o) 


(^t(x), y — X-|-^t(x)) otherwise. 


One can easily check that this equality implies Lemma For instance, for k >1, P(X = 1, y = 
k) = Q;a*(l — Q*)^ should be equal to 

P(X' = l,y' = A:) = = 1,X = 0,Y = k- 1) + P(6(x) = 1,X = l,y = A;) 

= p{l — a)a*{l — + paQ*{l — a*)^ 

= p{l — [(1 — a) + a(l — a*)] , 


which requires ([^. □ 

Theorem for Model 2 then follows from Lemma in the same way as before Theorem 
for Model 1 follows from Lemma [H □ 
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4 Law of Large Numbers in Hammersley’s processes 

In this section, we explain how Theorem on the stationary measure of the processes with 
sources and sinks implies Theorem 


Proof of Theorem\^ Trivial case. 


We first consider Model 1 in the case p > min{a/5, 6/a}. The upper bound is straightforward. 
The lower bound can be proven as follows. Let us consider for example the case p > ab~^. We 
have to prove: 


L 


lim inf 

n^oo 


( 1 ) 

{an,bn) 


> a. 


n 


( 12 ) 


One can build an increasing path as follows. The first point (l,yi) of the path is the lowest 
point of ^ having first coordinate equal to 1. The second point ( 2 , 1 / 2 ) of the path is the lowest 
point of ^ having first coordinate equal to 2 and second coordinate strictly larger than yi. The 
other points are defined in a similar fashion. As pb > a, a, study of this path easily provides 
(12). The case p > ba~^ can be handled in a similar way. 

The same strategy treats the case p > a/{a + b) in Model 2. 


Non trivial case 


We only need to compute lim since almost sure convergence follows from superad¬ 

ditivity. 

Model 1. Upper bound. We assume that p > mm{a/b,b/a}. For any a G (0,1), taking a* 
as in 0 and using ([^, ( |10| ) and we get 


^ between 1 and 6n}] 

= aa + ba* = aa + b —=: b, a). 

a+p{l-a) 


(13) 


The latter is minimized for 


a(o, 6) := 


VPVa-P 


a*{a, b) := 


1 ’ 

1 — P 

which both are in (0,1) if p < min{a/6, 6/a}. This yields 


VpVi-p 

1 ’ 

1 — p 


(14) 


Model 2. Upper bound. We assume that p > a/(a -|- 6). Using ( [TT| ) with a* = and 

that E[Geo(Q;*) — 1] = ^ — 1, we get 

{^^ks between 1 and 6n}] 

= aa + b{— -l) = aa + h^^^^—^. 

a* a — p 

We minimize the latter by taking 


a = p + y -p{l-p). 
a 


9 












Note that this choice is allowed since for p G (0, a /(o + h)) then 0 < p < a < 1. 

Model 1 and 2. Lower bound. 

In [9l Ho] the lower bound was obtained with a convexity argument on the scaling limit. 
Instead, we adapt a more probabilist argument due to Rolla-Sidoravicius-Surgailis-Vares (proof 
of Theorem 4.1 in [Hj) for the closely related model of last passage percolation with geometric 
weights. Here, we only give the details for Model 1 but the same argument applies for Model 2 
(there are only minor modifications due mainly to the fact that sinks have no more a Bernoulli 
distribution but a geometric distribution). 

Let us consider Model 1 on the rectangle [0, an] x [0, bn] with sinks and sources with optimal 
source intensity a = a{a,b) and a* = a*{a,b) (see ([l4])). For e G [0,1], denote by 
the length of the largest non-decreasing subsequence defined by 


= max{L; (ii,0),..., (ife,0), (4+i, jfc+i),..., (iL,jL) G I, 

0 < zi < • • • < ifc < ane < ik+i < ... < il < an and 0 < jk+i < • • ■ < jh < bn} . 

Here is an example where = 4. Note that because of the constraint is 

smaller than the number of Hammersley’s lines . 



Recall that has been defined in (13). 

Lemma 3 (Largest subsequence using some sources). There exist positive and non-decreasing 
functions f, g on (0,1], which depend on a,b, such that 


-/(e))) < exp(-ng(e)). 

( 1 ) 


Proof of Lemma 0 We write L^_ 6 „)(e) = Ii + /2 with 


Il := max{/c; (H,0),..., (4,0) G < • • • < 4 < one} 

h := max{/c; (4, ji),..., (4, jfc) G f,.,ane < 4 < • • • < 4 < an and 0 < ji < • • • < 4 < In} . 


The random variables 4 and 4 are independent, 4 is binomially distributed with parameters 
([anej, a(a, !))• random variable 4 has the law of I-‘^^an(i-£) bn)' dominated by a 


random variable 4 with law C 


(l),Q(a(l—£:),6),a*(a(l—e:),6) 
{an(l—e),hn) 


Using (13), we get: 



n 


pP) 

(an,i>n) 



< ^E[4+l3] 
n 

< aea(a, b) + ^a(l — e), b, a(a(l — e), 5)), 
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where, by convention, a{a',h') = 1 when p > a'/h' and a{a',b') = 0 when p > b'/a'. One can 
check that for p < niin{a/ 6 , b/a} and a', 6 ' > 0 such that a/b / a'/ 6 ' we have 

,b',a{a^b)) > (j)^^\a\b',a{a',b')). 


This is due to the fact that a(a,b) / a(a',b') and that there is a unique a G [0,1] which 
minimizes (/)(a', b', a). Thus, 


a£a(a, b) + 


^a(l — e),b, a(a(l — e), 6 )^ 


< aea{a,b) + ^a{l — e),b,a{a,b)^ 

= asa{a,b) + a{l — e)a{a,b) + ba*{a,b) 
= aa{a,b) + ba*{a,b) 

= (f>^^^a,b,a{a,b)). 


Denote by d{£) the difference between the right and left hand side of the inequality. The 
function d is positive on (0,1] and continuous. We set /(e) = ^ min{d((5), <5 G [e, 1]} which is 
non-decreasing and positive. Then, we have 

^ ^ ^ (^1 + ^3 > E(Ii + h) + n(d(e) - /(e))) 

< r{h + h> E(/i + h) + n/(e)). 


Recalling that Ii is binomially distributed and I 3 is the sum of two binomial random variables, 
function g is obtained by applying Hoeffding’s inequality to Ii and I 3 . 


□ 


We still consider Model 1 on the rectangle [0, an] x [0, bn] with sinks and sources with optimal 
source intensity a = a(a, 6 ) and a* = a*{a, b). Let 7 r„ be an optimal path for ■ If there 

are several optimal paths, we choose one of them in an arbitrary way. Define as the number 
of sources (z,/) G 7r„ with j = 0. 

Lemma 4 (Optimal paths do not take many sources). There exists a positive function h on 
(0,1] such that for and any 6 G (0,1] and for any n large enough 


¥ {Dn > and) < exp(—n/z((5)). 


In particular Dn/n converges to 0 in . 


Proof of Lemma By definition, if Dn > and then for some £ > d such that ean G N and e < 1 


• there are more than and sources in { 1 ,..., ean}; 


{an,on) ^ ’ {an,bn) 


'^(anb^ + if {sinks between 1 and bn}. 
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Therefore 


¥{Dn > anS} < P < n{aa - ^f{6)) 

+ P {sinks} < n{ba* - 

+ P bn)(^) ^ ^ “ fi^)) some e as above^ 

< exp{-ng{S)) + ^ ^ 

1>£>(5; eanGN 

< exp(-ng(S)) + ^ (^2i,bn) ^ 

1>£>(5; eanGN 

< exp(-ng(6)) + ^ exp(-n 5 r(e)) 

l>e>5; eanGN 

< exp(—ng(S)) + anexp(—ng(S)) 

for large re, where we used Lemma and where g is some positive function. This implies the 
first part of the lemma. The convergence of Dnjn to 0 in probability and in follows, since 
the sequence is bounded. □ 


b))-f{S))) 
b))- fie))) 


We now conclude the proof of the lower bound noticing that we have 

^(l),a a* ^ ^(1) 

(an,no) ^ — (an,on) 

where D'^ is the number of sinks (f,j) G tTu with f = 0. By the previous lemma we know that 
Dn/n tends to 0 in L^. By symmetry, the same result holds for D'^/n. Taking expectations of 
both sides in the previous inequality we get 

liminf-E[L[^)^^ J > aa + ha* = ff^\a,b,aia,b)). 
n-^+ca re (aiL,urb) 


□ 


Back to Ulam’s constant 

We observe that if we take p = 1/re, Theorem 1 for Model 1 suggests Ln^ ~ 2\/re, which is 
consistent with the asymptotics of Ulam’s problem, since ^ is then close, after renormalization, 
to a Poisson point process with intensity re. 

In fact, one can rigorously recover that c = 2 using our proof of Theorem 1. To do so, 
consider a Poisson point process H with intensity re in the unit square and denote by iin) the 
greatest number of points of H an increasing path can go through. To get a lower bound and an 
upper bound of £(re), we divide the square [0,1]^ in two different ways. 

First, we fix some A: > 1 and we divide it into small squares of length side 1/ik^/n). Say that 
^j{i) = 1 if at least one point of H is in the square with top-right corner [i/[k^/n), j/{ky/n)) 
and consider the quantity '■= ^-ssociated to the family iCji'i'))i,j<k^- is clear 

that 

((n) > (15) 

Denoting 

Pk = = 1) = P (Poiss(l/A:^) > l) = 1 — 
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Theorem implies 


^ky/n a.s. 2 /c-y/P^ 
-^/n n^+oo + 1 


Using (15) and letting k tend to infinity, we get c ;= lim£(n)/\/n > 2. 

To prove the upper bound, we divide now the square into small squares of length side 1/n^. 
Say that ij{i) = 1 if at least one point of H is in the square with top-right corner 
and consider the quantity := l |^4 ^ 4 ^ associated to the family The parameter 

of these Bernoulli random variables is now 


= P (Poiss(l/n^) > 1 ) = 1 — e 

With probability higher than 1 — n“^, all the columns and lines of width 1/n^ contain at most 
one point of H. On this event that we denote Fn-, coincides with i{n). We now use some 
intermediate results of the proof of the upper bound for a = 6 = 1. Inequality ( |10[ ) still holds 
despite the dependence of pn on n i.e. 


with = a* = ~ n 

"■ ^ l-Pn 


Using (13),we get 


< ■Y/E[Poiss(n) 2 ]/n 2 -|- -|- a*). 


where we used the Cauchy-Schwarz inequality. Dividing by n and letting n go to infinity, we 
obtain c < 2 . 


5 Stationary measures of Hammersley’s process on Z 

In this section, which can be read independently from the rest of the article, we study the 
analogous on the infinite line Z of the process {X^)t>o defined in ([^. This infinite Hammersley 
process was introduced by T.Seppalainen to prove the Law of Large Numbers for Problem 1. 

In order to generalize (W)t>o to a process taking its values in {0,1}^ we need a construction 
that does not rely on Hammersley lines. As before, (6(0)iezteN denote i.i.d. Bernoulli(p) 
random variables we say that there is a cross at time t located at x if ^tix) = 1. Informally, the 
infinite Hammersley process is defined as follows: 

At time t, if there is a particle at x and if the particle immediately on its left is at y, then the 

particle at x jumps at time t + 1 at the rightmost cross of f,t+i in interval {y,x) (if any, 

otherwise it stays still). 

For our purpose, we need a construction of (Xt) in which we let the crosses act one by one 
on configurations. We first draw Xq G {0,1}^ at random according to some distribution p such 
that ^—almost surely, for all i there is j < i such that A'o(j) = 1 . 

We now explain how to construct {Xt+i) from {Xt). A cross will act on a given configuration 
of particles in the following way: if a cross is located at x, then the leftmost particle of {Xt) in 
the interval [[x,-(-oo[[ (if any) moves to x; if there is no such particle then a particle is created 
at x. 

With this definition, we can now construct the value of (At+i) as a function of {Xt) and the 
crosses at time t -|- 1: We define Xt+i as the result of the successive actions on Xt of all crosses 
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at time t + 1 from the right to the left. Note that by definition, there only a finite number of 
crosses in that can modify Xt{i): those between Xt(i) and the first particle on its left. 

An example is drawn in this figure where the circles represent the particles and the crosses 
the locations where equals 1. Here, we have = (1, 0, 0,1, 0, 0,0,1,1,0,0,...) (as in the 
previous sections, time goes from bottom to top in our pictures): 



Theorem 3. The only extremal translation-invariant stationary measures of Hammersley’s pro¬ 
cess on Z are measures Ber(a)‘^^ for all a G (0,1]. 

This theorem is precisely the discrete counterpart of ([T], Lemma 7). The fact that Bernoulli 
product measures are stationary was proved by T.Seppalainen m Lemma 2.1). The reason for 
which we only focus on Model 1 on the infinite line is that the analogous of Model 2 is much 
more complicated to analyze. Indeed, the evolution of a particle in Model 2 depends on the 
whole configuration on its left. We don’t know if a similar statement to Theorem holds for 
Model 2. 

In order to prove that there are no other extremal measures we adapt the proof of Jockush- 
Propp-Shor ([5] p.17-20) for the discrete-time TASEP (see also [7] Ch.VIII). The main difference 
here with Jockush-Propp-Shor is our Lemma below that replaces their block argument. 

Proof of extremality. Let a G (0,1] and p be an extremal translation-invariant stationary mea¬ 
sure on {0,1}^ with marginals p{X{i) = 1) = a, we want to prove that p = Ber(Q;)®^. 

The proof is made of the following steps: 

1 . we introduce a measure w on { 0 , 1 }^ x { 0 , 1 }^ which is a minimal coupling between 
Ber(a)®^ and p; 

2. we prove (Lemma that some patterns that would decrease card{z;xj 7 ^ 7 /*}, have vr- 
probability zero; 

3. we conclude: p = Ber(Q!)®^. 

Let AI* be the set of translation-invariant measures on {0,1}^ x {0,1}^ that are a coupling 
of Ber(Q:)®^ and p. The set Af* is non-empty (it contains Ber(Q:)®^ ® p) and compact for the 
weak topology. 

We will prove that in Af* there is a coupling ((A(i))jgz, (h"(f))jez) such that X = Y a.s. 
Set 

D: ^ [0,1] 

TT ^ 7r(A(0)/y(0)). 

The function D is continuous on the compact set Af* and thus attains its minimum 5 > 0 at 
some coupling tt. If <5 = 0, then the theorem is proved: by translation-invariance X = Y 7 r-a.s. 

Let (XqjYq) ~ TT, we use the same random variables define the dynamics 

of X and Y. We denote by (Xt,Yt) the pair of configurations at time t, its joint distribution 
is denoted by vr^. By construction, for every t, is translation-invariant. For every A; > 1 we 
denote by 

Xk{Xt,Yt) 

the number of z G {1,..., A;} such that Xt{i) / Yt{i). 
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Let n > 2 and t > 0. We say that there is a n—forbidden pattern at location x G Z at time 
t if the configuration (Xt, Yt) between locations x and x + n — 1 is either 


X X + n — 1 

W = ( 1 0 ... 0 0 

y* = ( 0 0 ... 0 1 


) 

) 


or 


X X -\- n — 1 

W = ( 0 0 ... 0 1 ) 

y = ( 1 0 ... 0 0 ) 


Lemma 5. Forn > 2, let Ex^n be the event ^i{x) = 1, y(x+l) = y(x+2) = • • • = ^i(x+n—1) = 
0. Let j > 1, denote by F{n,j) the subset of locations of {0,n,2n, ■ ■ ■ ,{j — l)n} at which is 
located a forbidden pattern at time 0. Let 


A{n,j) = card{x G F{n,j), Ex^n occurs}. 


Then: 


^jn{Xi,Yi) < Ajn{Xo, Yq) + 1 — 2A{n,j). 


Proof of Lemma\^ We let the crosses at time 1 act one by one from the right to the left. We 
consider the impact of each cross action on the discrepancy Ajn between the two configurations. 
Let us note that, when a cross located at i acts on some configuration of particles, it changes 
the value of at most two sites of the configuration: i and the leftmost 1 of the configuration in 
the interval [[*,oo). Thus, we deduce the following facts: 


• The crosses strictly to the right of jn have no impact on the discrepancy. 


A cross located at a point x in {!,... ,jn} cannot increase the discrepancy. To check this 
fact, one can study all the cases. They are all shown below, up to symmetries between X 
and Y. 


/ X 

\ / 

X 



( 1 

) , ( 

1 

? • • 

• ? 1 ) 

V( 1 

) V( 

0 

0 •• 

•01)/ 




0 0 0 0 

0 1 ? ••• 




0 0 0 0 0 1 ) 

0 1 ? ••• ? 0 ), 


• For a forbidden pattern at some in, if Ein^n occurs, the cross at in decreases the discrepancy 
by 2. 

• In the worst case, the right-most cross which is strictly to the left of 1 increases the 
discrepancy by 1. 

• The crosses to the left of the previous one have no impact on the discrepancy. 

The result follows. □ 


Lemma 6 (Forbidden patterns). Let n >2. The ttq probability of a n-forbidden pattern at any 
location is 0. 
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Proof of Lemma^ Let j > 1. By Lemmawe have: 


^jn{Xi,Yi) < Ajn{Xo, Yq) + 1 — 2A{n,j). 

Note that the probability that, at some location in, there is a n-forbidden pattern and that 
occurs is Cp( 1 — where C is the probability of a n-forbidden pattern at a given location. 

We aim to prove C = 0. Taking expectation in the previous display we get 

jnDfe) < jnD{-KQ) + 1- 2jCp{l 

Dividing by nj and letting j —>■ oo we get: 

2Cp{1-p)^-^ 


D{Ki) < D{7Lo) - 


n 


But the minimality of tt yields D{ni) > D{t^q) and we get C = 0- 


□ 


Lemma 7 (There is a forbidden pattern somewhere). Assume that <5 > 0, there exists n such 
that 


> 0 



1 

n \ 


1 

n \ 

x = 

( 1 ? . 

.. ? 0 ) 

>0 or TT X = 

( 0 ? . 

.. ? 1 ) 

V y = 

( 0 ? . 

.. ? 1 ) y 

' V ^ = 

( 1 ? . 

.. ? 0 ) / 


Proof. Let {X, T) ~ tt, 

5 = 7r(x(o) / y(o)) = 7r(x(o) = 1 , y(o) = 0 ) + 7r(w(o) = 0 , y(o) = 1 ). 

The two terms on the right side are equal: 


( 0 ) = 1 , y(o) = 0 ) = ^(x(o) = i) - ^(w(o) = i, y(o) = i) 

= 7r(y (0) = 1) — 7r(y (0) = 1, y (O) = l) {X, Y have same marginals) 

= ^(x(o) = o,y(o) = i). 

Then 0 < 5/2 = ;^(X(0) = 1, y(0) = 0). Assume that the lemma is false (for every n), we have 
7r(A U B) = 1, where 


A = {{x, y), Xi < Vi for all i G Z} , 

B = {{x,y),Xi> yi for all i G Z} . 

As seen above, 5 > 0 implies then that 7r(A),ff(i?\A) > 0. Besides, it is easy to check that A 
and B\A are preserved by the dynamics. Hence, considering the restriction of on these two 
subsets, := {va, Pa) and tIb\a ■— {^b\Aj Pb\a)^ see that va (resp. pa) and u^^a (resp. 
Pb\a ) define two translation-invariant stationary measures such that 

Ber(a)®^ = 7r{A)nA + 7 l{B\A)ub\a- 
P = 1L{X)pa + 7l{B\A)pb\A- 

The extremality of Ber(a)®^ and p implies that Ber(Q;)®^ = i^A = t^b\a and p = pA = Pb\a- 
But, by definition of A and B, va ^ PA and Pb\a ^ ^B\A- Thus, necessarily, p = Ber(Q;)®^. □ 
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Now we can obtain our contradiction. Indeed, there exists an integer n such that 



1 

n \ 


1 

n \ 


( 1 ? . 

.. ? 0 ) 

\ > 0 or vr X = 

( 0 ? . 

.. ? 1 ) 

V = 

( 0 ? . 

..?!)) 

V ^ = 

( 1 ? . 

.. ? 0 ) / 


Without loss of generality we do the first case. We can assume that X, Y coincide between 
2 and n — 1, otherwise we can decrease n. Let 1 < A: < n be the leftmost position at which 



With positive probability we can turn this (j^) into (J): 


/ 


X — - 





f 





\ 



1 

k 

n 


t+i 



1 

k 

n 



w = 

(10.. 

. 0 1 

? 0 

) 



Xt+i = 

( ? 0 .. 

.01? 

0 

) 

V 

Yt = 

(00.. 

. 0 1 

? 1 

) y 


V 

rt+i = 

( ? 0 .. 

.00? 

1 

) / 


We repeat this process between positions k and n and finally for some t,i,j we have 

/ ^ J \ 

vr Xt= ( 1 0 ... 0 0 ) >0, 

V = ( 0 0 ... 0 1 ) / 

which contradicts Lemma [H □ 


A 


Dynamic of the underlying interacting 


particle system 


In this appendix we explicit the dynamic of the underlying interacting particle system in Models 
1 and 2. Consider the processes {Xl)t>o defined by Q (without sources and sinks). One can 
explicitely construct the dynamic of these processes in term of the family of i.i.d. random 
variables {^t{x),x G [[l,n]],t G N} with law Bernoulli(p) (see Q). To do so, let say that there is 
a cross at time t located at x if ^tix) = 1. A cross will act on a given configuration of particles 
in the following way : if a cross is located at x, then the left-most particle in the interval [x, nj 
(if any) moves to x ; if there is no such particle then a particle is created at x. 

With this definition, we can now construct the value of as a function of XI and the 
crosses at time t + 1: 


Model 1. We define as the result of the successive actions on X^ of all crosses at time 

t -|- 1 from the right to the left. 

An example is drawn in this figure where the circles represent the particles and the crosses 
the locations where ^t+i equals 1. Here, we have = (1, 0, 0,1, 0, 0,0,1,1,0,0) (as it is usual 
in the literature on Hammersley’s processes, time goes from bottom to top in our pictures): 





C) 

o 






n 


Let us note that there is an alternative way to construct the model 1 : At time t, if there is 
a particle at x and if the particle immediately on its left is at y, then the particle at x jumps 
at time t -|- 1 at the rightmost cross of in interval (y, x) (if any, otherwise it stays still). If 
there are some points of ft+i at the right of the rightmost particle of Xf then a new particle 
appears at the leftmost point among them. 

Model 2. We define as the result of the successive actions on A| of all crosses at time 

t + 1 from the left to the right. 
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We let the reader convince themselves that this construction coincides with the definition 
given in Q. 

Although the definition of both models seems, at first glance, very close, the nature of the 
two processes is in fact quite different. Indeed, in the hrst model, to find the location at time 
t + 1 of a particle located at time t at x, one just need to know the location y of the particle 
immediately on his left at time t and the position of the crosses of in the interval {y, x). In 
particular, with no difficulty, one can define a process with similar transitions on the whole line 
Z as it is done in [9]. For Model 2, the dependances are more intricate. Indeed, to determine 
the location at time t + 1 of a particle located at time t at x, one need to know the whole 
configuration of Xf and ^t+i on the interval [[l,x]]. In particular, the definition on the whole 
line Z of a similar process is more delicate and requires a condition between the density of crosses 
and particles (see [T0]b 

Finally, let us note that the construction can be extended to the models with sources and 
sinks. The set of sources simply give the initial configuration of the particle system. The sinks 
acts as the crosses. In Model 1 sinks act after the other crosses while in Model 2 sinks act before 
the other crosses. 
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